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A GOODNESS-OF-FIT  TEST  FOR  A CLASS  OF  NONHOMOGENEOUS 
POISSON  PROCESSES 


1.  INTRODUCTION 

The  problem  addressed  here  is  that  of  testing  the  goodness- 
of-fit  of  certain  integer  valued  stochastic  processes.  The  hypothesis 
to  be  tested  is  that  a set  of  waiting  times  are  from  a nonhomogeneous 
Poisson  process  which  is  a member  of  the  class  of  processes  with  mean 
value  function  of  the  form 


M(t)  = Xt®  t>0.  (1.1) 

The  occurrences  of  failures  in  military  systems  undergoing  refinement 
in  a test-fix-test-fix  development  process  and  in  complex  systems  which 
are  repaired  upon  failure  are  known  to  follow  such  processes.  Crow  (1) 
has  discussed  this  application  and  shown  how  to  use  the  Cramdr-von  Mises 
statistic  to  test  the  goodness-of-fit  hypothesis.  A table  of  the  small 
sample  distribution  of  the  Cramdr-von  Mises  statistic  for  the  case  in  which 
an  exponentially-appearing  parameter  is  estimated  also  appeared  in  (1). 

Some  percentiles  of  that  table  are  inaccurate  because  of  sample  sizes  used 
in  the  simulations  which  generated  them.  The  objectives  of  the  present 
work  are  to  provide  a more  accurate  table  and  to  verify  convergence  to 
the  asymptotic  distribution.  In  addition  various  proofs  were  omitted 
from  the  previous  work.  Some  of  those  proofs  are  included  here  for 
completeness . 

The  procedures  presented  here  can  also  be  used  to  test  the 
hypothesis  that  a random  sample  from  a continuous  distribution  has  a 
cumulative  distribution  function  F(X;0)  of  the  form 

F(X;0)  = (R(X))®  (1.2) 

for  some  positive  value  of  the  parameter  0 and  for  a specified  cumulative 
distribution  function  R(X).  Tables  of  the  small  sample  distribution  of 
the  Cramdr-von  Mises  test  are  provided  to  implement  the  test. 

A 

2.  THE  CRAMER-VON  MISES  TEST  IN  THE  PARAMETRIC  CASE 

Cramdr  (2),  von  Mises  (3),  and  Smirnov  (4)  developed  a test 
of  the  hypothesis  that  a random  sample  X^  X2,.  . . ,Xj^  from  a continuous 

distribution  G(X)  is  drawn  from  the  completely  specified  distribution 
function  F(X).  If  the  cumulative  distribution  function  F(X)  contains  a 
parameter  0,  then  the  hypothesis  to  be  tested  is  that  G(X)  = F(X;0  ^ 

o'* 

for  some  specified  0^.  The  statistic  employed  in  the  test  is  given  by 

[F^(X)-F(X;0^)]2  dF(X;0^). 

ft#  mOO 
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(2.1) 


where  F (X)  is  the  empirical  distribution  function.  The  empirical 
distribution  function  is  defined  as  if  k of  the  Xj^  are  less 

than  X.  It  can  be  shovm  that  the  statistic  above  can  be  written  as 
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where  X.  is  the  j 'th  order  statistic  of  the  random  sample, 

3 2 

that  G(X)  = F(X;0  ) is  rejected  if  W is  extraordinarily  large 


The  hypothesis 
Smirnov 


derived  the  limiting  distribution  for  as  the  sample  size  n becomes 
large.  Anderson  and  Darling  (5)  tabulated  this  distribution. 


Following  a suggestion  of  Cramdr,  Darling  (6)  extended  this  test 
to  the  case  in  which  the  parameter  0 is  estimated  by  a statistic 

0_  calculated  from  the  data.  This  test  assumes  that  there  exists  a 
n 

nondegenerate  interval  I on  the  real  axis  such  th'^t  for  every  0 contained 
in  the  interior  of  I,  F(X,'0j  is  a cumulative  distribution  function.  The 
hypothesis  tested  is  that  the  cdf  from  which  the  sample  is  drawn  is  a 
member  of  the  parametric  family  F(X;0)  for  some  unknown  value  0 in  I. 


For  the  parametric  case  the  test  statistic  is 
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F(X.;0^) 
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in  which  the  estimate  0^  is  substituted  for  the  unknown  parameter.  Darling  (6) 
investigated  the  limiting  distribution  of  and  found  that  there  are 
essentially  two  distinct  cases.  The  first  case  is  that  of  a superefficient 
estimator  0,  such  that 


lim  nE  =0.  (2.4) 

n-*- » ' 3 

where  0 is  the  true  value  of  the  unknown  parameter.  In  this  case  the 

2 2 

asymptotic  distribution  of  C is  the  same  as  that  of  W . In  the  second 

n 2 ^ 

and  more  general  case  the  limiting  distribution  of  C is  different 
2 

from  that  of  W . Thus  it  is  inappropriate  to  use  the  tabled  critical 


values  for  '.v'“  to  test  the  gocdness-of-fit  h^'pothesis  in  this  case.  If 

certain  regular! tv  conditions  are  satisfied  and  id  -d''  is 
' ■ n ' 

asyTTipotically  nortial  '.vith  xean  tero  and  finite  nontero  variance,  then 

a . 1 ' 

the  limiting  distribution  of  C“  is  the  same  as  that  of  r * Y'ft)  dt 

where  Y(t)  is  a Gaussian  stochastic  process  with  cero  mean  and  a 
covariance  kernel  which  depends  upon  the  unknown  true  distribution 

7 

G(X).  Thus  the  C"  test  is  not  distribution  free;  however,  in  important 
special  cases  the  aistribution  depends  only  upon  the  form  of  the  family 
f(X,9)  and  not  upon  the  true  value  9^.  Fortunately,  for  the  case  of  an 

e.xponential  parameter  which  is  of  concern  in  this  report,  the  test  has 
this  property  of  being  parameter  free. 

3.  NON'HOMOGENOUS  POISSON  PROCESSES 

Parcen  (7)  defines  the  concept  of  a nonhomogeneous  Poisson 
process.  This  is  an  integer  valued  process  which  has  independent 
increments  and  unit  jumps  and  which  consists  of  the  random  variables 
NTt),  the  number  of  events  occurring  in  the  interval  I0,t].  It  follow-'? 
from  this  definition  that  for  any  t>0,  the  random  variable  X'Ct)  has 
the  Poisson  distribution  with  an  e.xpected  value  given  by  some  function 
MCt).  The  mean  value  function  M(t)  is  nondecreasing  and  is  usually 
assumed  to  be  continuous  and  differentiable.  The  derivation  of  the 
mean  value  function  is  denoted  by 


vCt)=^MCt)  (3.1) 

and  is  called  the  intensity  function  of  the  process.  For  a small  in- 
crement of  time  h the  quantity  v(t)  h is  approximately  equal  to  the 
probability  of  the  occurrence  of  an  event  in  the  interval  (t,t+h). 

In  the  special  case  for  which  the  mean  value  function  is 
directly  proportional  to  the  observation  time  t,  that  is, 

MCt)  = Xt  (3.2) 

the  process  is  the  ordinary  homogeneous  Poisson  process.  Such  a pro- 
cess has  a constant  intensity  function. 

In  the  general  case  of  the  nonhomogeneous  Poisson  process 
use  of  the  characteristic  function  along  with  the  property  of  indepen- 
dent increments  reveals  that  for  any  interval  (a,b)  the  number  of  events 
occurring  in  the  interval  has  a Poisson  distribution.  The  characteristic 
function  of  N(t)  is 


'^‘NCt)^^^  = explMCt)  {e^'^  -1}] 


(3.3) 


’I 


Since  the  numbers  of  events  occurring  in  nonoverlapping 
independent  we  have 


"N(b)-NCa) 


(u)  _ ^NCb)  _ 


= exp  ( [MCb) -M(a)  ] [e^^-1 


NCa) 


intervals  are 


] . C3.4) 


Thus  the  random  variable  N(b)-NCa)  has  the  Poisson  distribution  with 
expected  value  M(b)-M(a). 

A pair  of  well  known  theorems  describe  the  relationship 
between  the  order  statistics  of  a random  sample  from  a uniform 
distribution  and  the  waiting  times  to  the  occurrence  of  events  in  a 
homogeneous  Poisson  process.  The  proof  of  the  first  is  given  by  Parzen. 
The  statements  of  both  theorems  and  the  proofs  are  presented  here  to 
clarify  this  relationship. 

Theorem  1.  Let  { N(t),  t^O}  be  a Poisson  process  with  constant  intensity 
X.  Under  the  condition  that  N(T}=n,  the  times  t^,t2,...t^  in  the 

interval  (0,T]  at  which  events  occur  have  the  joint  probability  density 

function 


r 

t 


1 


I 

f 

t 

t 


fCti.t2.---.tn 


N(T>n)  = 


(3.5) 


in  which  0 < t^  < t.,  joint  conditional 

density  is  the  same  as  that  of  the  order  statistics  of  a random  sample 
size  n drawn  from  a uniform  distribution  defined  on  the  interval  (0,T]. 


Proof;  The  probability  that  exactly  n events  occur  in  the  interval 
(O.T]  is 

P[NCT)=n]  = . C3.6) 

nl 


Let  {t.}  be  a set  of  times  such  that  o <t,  <t_<* • • <t  < T and  choose  a 
1 1 n 

set  of  increments  fli^,h2,  • • • ,h^}  such  that  the  intervals  ■••hj], 

[t2,t2+h2] , • • • [tn,tn+hn]  are  nonoverlapping.  Since  the  increments  are 

nonoverlapping  the  number  of  occurrences  in  each  interval  and  the 
number  occuring  elsewhere  in  CO>T]  are  independent.  Thus  the  probability 
that  exactly  one  event  occurs  in  each  interval  [t.,t.+h.]  for  ial,2...,n 
and  no  events  occur  elsewhere  in  (0,T]  is  i i i 


1 


4 


n 

e.\p[-  \ (T - 'i'  h. ) ] 

i = l * 

The  conditional  .probability  of  this  event  given  that  \'(t)  = n is 
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i=l  " 
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3y  .iliowing  the  h.  to  decrease  to  Lnfin itesir.al  increments,  it  follows 
that  the  conditional  density  of  t^  ,t.,, . . . ,t^  given  that  N'(t)=N  is  as 
stated  in  the  theorem. 


Iheorem  d . Let  !;.\;_t  j .t^^O;  be  a foisson  process  with  constant  intensity 

Let  t be  the  time  of  occurrence  of  tiie  nth  e\'ent . under  the 
n 

condition  that  t =T,  the  n-1  waiting  times  t,,t,,...,t  , have  the  joint 

n =■  1 2 n-1  ^ 

probability  density  function 


fCt^.t^, . . . ,t 


n-1 


t =T)  = 
n .^n-1 


C3.9) 


N’ote  that  this  point  conditional  density  is  the  same  as  that  of  the 
order  statistics  of  a random  sample  of  sice  n-1  drawn  from  a uniform 
distribution  defined  on  the  interval  C0,T]. 

Proof;  Define  t =0  and  let  s.=t.-t.  , be  the  i-th  interarrival  time 
0 11  1-1 

of  the  Poisson  process.  It  is  kno\m  that  the  interarrival  times  are 
independently  identically  distributed  random  variables  with  the  e.x- 
ponential  distribution  with  mean  1/X.  Thus  the  joint  density  of  the 
first  n interarrival  times  is 

n -Xs.  n 

f(s  ,s  ) = IT  Xe  ^ = x’’  exp[-X  ) s ] (5.10) 

1=1  1=1 

for  which  0 <s^ , s.^, . , . , s^  1’^"'  transformation  mapping 

(tj;  ,t,, . . . ,t^  into  (s^.s.,,...,  ^ has  a Jacobian  identically  equal 

to  unity.  Therefore  the  joint  density  of  the  waiting  times  is 


gCt, ,t 


l’"2’' 


- 


x"e 


-Xt 


5 


(3.11) 


with  0<t,  . . . <t  . 

1 _ n 


t 15  tr.-?  convol’iticr’.  or  n rar.^oi  '.'iri'i’''!  r""' 

n 

sane  exponent iai  distribution,  it  has  the  gatiita  density 


g [t  ) = 
*n  n 


■■’^nCXtn)'^-' 

Cn-1)! 


(3.12) 


It  follows  taat  the  conditional  density  of  t,,...,t 

t =T  is  as  in  the  statement  of  the  theorem, 
n 


n-1  given  tr.at 


A nonhomogeneous  process  with  a continuous  mean  value  function 
.'li.t)  can  be  transformed  into  a homogeneous  .-'oisscn  process,  because 

'•l(t)  is  continuous  and  nondecreasing  its  inverse  function  ^ (X)  car.  be 
defined  for  all  X^O  as  the  minimum  value  of  t such  that  M(t)^X.  Define 
the  stochastic  process  '’xfX)  ,X>  O'- 


X(X)  = \'Cm'\'X)) 


(3 . 13) 


This  is  a Poisson  process  with  an  intensity  function  identically  equal 
to  one.  This  transformation  will  be  used  in  the  proof  of  the  two  sub- 
sequent theorems. 

Theorem  3.  Let  '>;(t),t^0)-  be  a Poisson  process  with  continuous  mean 

value  function  MCt).  Under  the  condition  that  N(T)  = n,  the  n waiting 

times  tT,t.,,...,t  in  the  interval  (O.T]  at  which  events  occur  are 
1 ^ n 

random  variables  having  the  same  distribution  as  the  order  statistics 
of  a random  sample  of  size  n from  the  probability  density 


hfXX  * 

“ MTfT 


0<X<  T 


C3.14) 


where  v(X)  = dM(X)/dx  is  the  intensity  function  of  the  Poisson  process, 
That  is  , 


gCtj.t2,...,t^  NCT)=n) 


with  0<t,  , <t.<.  . . <t  <T. 
— 1 2 n— 


n ^ 

[MCT)]  i=l 


(3.15) 


L 


Bf  /v- 


Proof;  Define  the  inverse  of  the  mean  value  function  and  the 
stochastic  process  {KCX),x>0'  as  in  the  preceding  paragraph.  It 
that  the  ouantities  MCt.)  = X.  for  i = are  the  first  N' 

1 i 


foilous 

wait- 


ing tires  fro.u  a Poisson  proces;  ..ith  intensity  identioally 
unity  and  occurring  in  the  interval  CO/'iT)].  3y  theorem  1 
conditional  densitv  of  X,,...,X  given  that  N[MfT)]  = n is 

' i .X 


^ciual  to 
t h ^ i 0 i n t 


fix. ,x. 


^ MfT'i  >=• 


The  absolute  value  of  the  Jacobian  of  the  transformation  which  tiaps 

(t,  ,t,, . . . ,t  ) into  (X, ,X, X ) is  given  by 

1 - n i - i4 

I I 


i = l it. 


, which  is 


J|  = P vCt  ) 
i = l ^ 


(5.1-) 


It  follows  that  the  conditional  density  g (t^  , t , . . . , t^  jX (t  j =n)  is  as 

given  in  the  theorem.  Note  that  the  distribution  function  corresponding 
to  the  density  h(X5  is  given  by 


H(X) 


M(X1 

M(T) 


0<X<T 
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Theorem  4 . Let  {N(t),t>0}  be  a Poisson  process  with  continuous  mean 
value  function  MCt).  Under  the  condition  that  the  w'aiting  time  to  the 
nth  event,  t^,  is  equal  to  T,  the  n-1  waiting  times  t j^<t,<.  . . <t^  in 

the  interval  (0,T)  at  which  events  occur  are  random  variables  having 
the  same  distribution  as  the  order  statistics  of  a random  sample  of 
size  n-1  from  the  probability  density 

h(X)  = ^ 0<X<T  C3.19) 


where  v(X} 
That  is, 


dMC.x) 

dx 


is 


the  intensity  function  of  the  Poisson  process. 
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wit!'.  0<t,<T,". 


•'  i extent  that  tl'.eorem 
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conternir.';  the  nean  value  fv.netion  of  a voi 
oarar.eter  v.  as  ‘ollcws 
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:or  tae  eor.a.tier.  .'h  • =a 


n-1  for  the  condition  t =T 

n 


K''  * - ~ J 


in  which  n is  the  nurtber  o^  events  obser''ed  in  che  interval  rO,T]  of  a 

stochastic  nrocess  (XCt' , t>0' . Consider  t'ne  hvnothesis  H r-Ctj  is  a 

— • - 0 

Poisson  process  with  continuous  mean  value  function  'l(t) . If  is 

true  the  waiting  times  t^,t^,...,t  have  the  same  distribution  as  the 

order  statistics  of  a random  sample  of  sice  m from  the  distribution 
H(;X}=M(X] /MCT)  . Define  H^(X)  to  be  >I(X)/N(T).  Then  the  statistic 


'"O  [H_^CX)-HCX)]"dHCX) 


(5.25) 


has  the  same  distribution  as  the  Cram^r-von  Mises  statistic  for  a 
sample  of  sice  m from  HfX) . The  statistic  can  be  written  in  the  form 


, 1 ™ ''(t.)  2.  , , 

w-  = ^ - y r L . 1 “ 

m 12m  .-/'-'(T)  ' ’ 

1 = 1 2m 


C5.24'i 


AVA!L/;B; . 


which  is  more  suitable  for  computation. 

If  the  mean  value  function  contains  an  unknown  parameter  0, 
then  it  is  desirable  to  estimate  the  parameter  from  the  data  by  calculat- 
ing a statistic  0^^^.  If  the  estimator  0^  satisfies  the  properties 
listed  by  Darling,  then  the  statistic 
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12m 
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j = l 


M(T.Om) 


(3.25) 


may  be  used  to  test  the  hypothesis  : {N(t)t^0}  is  a Poisson  process 

with  mean  value  function  MCT;0)  for  some  0.  The  test  is  truly  usable 
2 

if  the  statistic  is  parameter- free . 

4.  THE  RELIABILITY  GROWTH  PROCESS 

Crow  (8)  has  shown  that  the  improvement  in  reliability  of  a 
complex  system  undergoing  development  in  a test-fix-test-fix  environment 
can  be  modeled  by  a certain  family  of  nonhomogeneous  Poisson  processes. 
Crow  (1)  uses  this  same  family  of  processes  to  represent  the  occurrences 
of  failures  in  complex  repairable  systems.  For  this  class  of  processes 
the  mean  value  function  is  of  the  form 

M(t)  = \t^  X>0;6>0;t>0  (4.1) 


in  which  X>0  can  be  interpreted  as  a scale  parameter  and  6>0  as  a shape 
parameter.  The  corresponding  intensity  function  is 

v(t)  = X6t^'^  t>0  . (4.2) 


This  family  includes  the  homogeneous  Poisson  processes  as  the  snecial 
case  in  which  6 equals  unity. 

The  results  of  the  preceding  section  can  be  used  to  derive  a 
goodness  of  fit  test  for  this  class  of  processes.  With  the  index  m 
defined  as  in  Section  3,  it  follows  from  theorems  3 and  4 that  with 
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application  of  the  appropriate  condition  the  random  variables  >^2 > • • • 

have  the  same  distribution  as  the  order  statistics  of  a random  sample 
for  the  cumulative  distribution  function 

H(X)  OlX£r  . (4.3) 

With  the  appropriate  choice  of  a estimator  6^  the  statistic 


C 


2 

m 


1 

12m 


2 


(4.4) 


can  be  used  to  test  the  hypothesis  that  the  observations  (t.)  are  from 

J ft 

Poisson  process  with  mean  value  function  of  the  form  M(t)  = At  • Darling 
has  shown  that  the  distribution  of  this  test  statistic  is  independent 
of  the  true  value  of  the  parameter  6.  In  fact  it  is  distribution  free 
over  the  class  of  distributions  such  that 


F(X;6)  = [R(X)j  ® B>0  (4.5) 

for  some  cumulative  distribution  function  R(X). 

As  Crow  has  shown  conditional  maximum  likelihood  estimates  of 

B can  be  derived  from  equation  (3.15)  for  the  case  N(T)  = n and  from 

equation  (3.20)  for  the  case  t = T.  These  conditional  maximum 

n 

likelihood  estimates  are  given  by 


m 

_ 

mint  - i ^i 
i=l 


(4.6) 


in  which  m = n for  conditioning  on  N(T)  = n and  m = n-1  for  conditioning 


This  is  a biased  estimator  for  6 with  expected  value 


on  t » T. 
n 


E(?i)  - 6 

m*  A 


(4.7) 


An  unbiased  estimate  can  thus  be  provided  by 


m-1 

m 

mini  - '^i 

i<=l 


(4.8) 


It  can  be  shown  that  the  estimator  8 given  in  equation  (4,8) 

2 

should  be  used  in  equation  (4.4)  to  calculate  the  statistic  C . This  must 

2 ro 

be  done  in  order  to  satisfy  the  conditions  needed  for  Cj^  to  have  the 

limiting  distribution  described  by  Darling.  In  particular,  the 
condition 


lim  mE(6-8)  = 0 (4.9) 

is  satisfied  since  F is  unbiased.  This  condition  is  not  met  by  the 
estimator  in  equation  (4.6). 


5.  DISTRIBUTION  OF  THE  STATISTIC 

m 

2 

In  order  to  use  the  statistic  to  test  the  goodness  of  fit 

hypothesis  it  is  necessary  to  establish  a table  of  critical  values  for 
selected  significance  levels.  The  small  sample  distribution  of  the 
2 

statistic  given  in  equation  (4.4)  is  not  analytically  tractable. 

Moreover,  the  limiting  distribution  has  only  been  defined  in  terms  of 

2 

its  characteristic  function.  The  distribution  of  has  been  determined 

through  Monte  Carlo  simulation  for  values  of  m from  2 to  20  and  for 
m equal  to  30,  60,  and  100. 

The  Monte  Carlo  simulation  consists  of  repeated  samples  of 
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size  m from  the  uniform  distribution  on  the  interval  (0,1)  and 
computation  of 


6 = 


m-1 


nt 

-I  In  U 

3 = 1 ^ 


(5.1) 


and 


< ■ ik  • I 

3 = 1 


2j-l 

2m 


(5.2) 


in  which  Is  the  random  sample  and  (uj)  is  the  corresponding  set  of 

order  statistics.  For  each  value  of  the  index  m there  are  150,000 
replications  of  this  sampling. 

2 

Selected  percentage  points  of  the  distribution  of  are 

presented  in  Table  1.  The  1-a  percentile  of  this  distribution  is  to  be 
used  for  a goodness  of  fit  test  with  level  of  significance  a.  The 
accuracy  of  these  percentage  points  can  be  determined  by  using  the  fact 
that  any  percentile  of  a random  sample  is  asymptotically  normal.  Each 
sample  of  150,000  actually  consists  of  ten  independent  samples  of  size 

15,000.  The  sample  variance  of  estimate  C^.  of  the  p-th  percentile  is 

used  to  estimate  the  precision  of  the  p-th  percentile  of  the  combined 
sample.  Table  2 contains  interval  estimates  of  the  percentiles  of  the 
2 

distribution  of  C_,  with  a confidence  coefficient  of  .90. 
m 


The  sample  moments  from  the  simulation  can  be  used  to  determine 

2 

how  rapidly  the  distribution  of  is  converging  to  the  limiting 

distribution  . Darling  provided  a means  for  calculating  the  moments 
of  the  limiting  distribution.  The  mean,  variance,  and  third  central 

2 

moment,  of  the  sampling  distribution  of  and  of  the  limit 

distribution  are  given  in  Table  3.  The  sample  moments  indicate  the 
distribution  for  m equal  to  100  matches  the  limiting  distribution  quite 
closely.  The  mean,  variance,  and  third  central  moment  from  the  simula- 
tion for  m equal  to  100  are  each  within  one  percent  of  the  respective 
true  moment  of  the  limiting  distribution.  Hence  the  percentiles  appearing 
in  ’.’able  1 for  m equal  to  100  can  be  used  for  larger  values  of  m.  Figure 
1 contains  plots  of  the  empirically  obtained  density  function  for  m equal 
to  5 and  m equal  to  100. 


12 


6. 


CONCLUSION 


The  percentiles  in  Table  1 provide  a set  of  critical  values 
for  the  Cramdr-von  Mises  goodness  of  fit  statistic  for  the  case  in  which 
an  exponential  parameter  is  estimated.  This  table  can  be  used  to  test 
the  hypothesis  that  a random  sample  of  size  m comes  from  a parametric 
family  in  the  class  of  distribution  of  the  form 


F(X;0)  » lR(X)l  0>O  (6.1) 


in  which  RCX)  is  some  cumulative  distribution  function.  The  parameter 
0 is  to  be  estimated  from  the  data  by  an  appropriate  statistic. 

The  distribution  of  the  Cramdr-von  Mises  statistic  can  also 
be  used  to  test  hypotheses  on  the  goodness  of  fit  for  certain  stochastic 
processes.  In  particular,  the  hypothesis  that  a stochastic  process  is 
a member  of  the  family  of  nonhomogeneous  Poisson  processes  with  mean 
value  function  of  the  form 


MCt)  = Xt® 


(6.2) 


can  be  tested  through  use  of  the  statistic. 


TABLE  1 

PERCENTILES  FOR  THE  DISTRIBUTION  OF 


O 

00 

.85 

.90 

.95 

.99 



2 

.138 

.149 

.162 

.175 

.186 

3 

.121 

.135 

.154 

.184 

.23 

4 

.121 

.134 

.155 

.191 

.28 

5 

.121 

.137 

.160 

.199 

.30 

6 

.123 

.139 

.162 

.204 

.31 

7 

.124 

.140 

.165 

.208 

.32 

8 

.124 

.141 

.165 

.210 

.32 

9 

.125 

.142 

.167 

.212 

.32 

10 

.125 

.142 

.167 

.212 

.32 

11 

.126 

.143 

.169 

.214 

.32 

12 

.126 

.144 

.169 

.214 

.32 

13 

.126 

.144 

.169 

.214 

.33 

14 

.126 

.144 

.169 

.214 

.33 

15 

.126 

.144 

.169 

.215 

.33 

16 

.127 

.145 

.171 

.216 

.33 

17 

.127 

.145 

.171 

.2n 

.33 

18 

.127 

.146 

.171 

.217 

.33 

19 

.127 

.146 

.171 

.217 

.33 

20 

.128 

.146 

.172 

.217 

.33 

30 

.128 

.146 

.172 

.218 

.33 

60 

.128 

.147 

.173 

.220 

.33 

100 

.129 

.147 

.173 

.220 

.34 
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TABLE  2 


INTERVAL  ESTIMATES  OF  PERCENTILES  OF  THE  DISTRIBUTION 
90%  CONFIDENCE  COEFFICIENT 


.80 

.85 

.90 

.95 

.99 

M X 

2 

.1372-. 1378 

.1485-. 1493 

.1613-. 1618 

.1748-. 1753 

.1863-. 1866 

3 

.1203-. 1212 

.1346-. 1356 

.1541-. 1548 

.1828-. 1843 

.2293-. 2320 

4 

.1204-. 1211 

.1337-. 1346 

.1539-. 1552 

.1901-. 1924 

.2761-. 2819 

5 

.1202-. 1212 

.1359-. 1372 

.1588-. 1605 

.1984-. 2006 

.2923-. 2978 

6 

.1220-. 1232 

.1384-. 1395 

.1609-. 1633 

.2018-. 2054 

.3029-. 3115 

7 

.1233-. 1242 

.1396-. 1406 

.1639-. 1652 

.2062-. 2096 

.3128-. 3202 

8 

.1233-. 1245 

.1400-. 1414 

.1644-. 1663 

.2084-. 2106 

.3143-. 3238 

9 

.1245-. 1253 

.1419-. 1428 

.1667-. 1683 

.2109-. 2130 

.3191-. 3261 

10 

.1243-. 1257 

.1416-. 1429 

. 1663- . 1684 

.2108-. 2136 

. 3208- . 3264 

11 

.1253-. 1263 

.1429-. 1439 

.1681-. 1691 

.2128-. 2160 

.3249-. 3302 

12 

.1260-. 1268 

.1431-. 1441 

.1679-. 1697 

.2124-. 2151 

.3215-. 3276 

13 

.1263-. 1272 

.1440-. 1451 

.1695-. 1711 

.2140-. 2164 

.3230-. 3312 

14 

.1259-. 1269 

.1435-. 1445 

.1686-. 1697 

.2130-. 2154 

.3261-. 3322 

15 

.1256-. 1270 

.1433-. 1451 

.1680-. 1707 

.2138-. 2164 

.3257-. 3301 

16 

.1266-. 1278 

.1440-. 1454 

.1701-. 1714 

.2146-. 2174 

.3216-. 3291 

17 

.1268-. 1281 

.1446-. 1455 

.1697-. 1713 

.2152-. 2181 

.3288-. 3351 

18 

.1271-. 1278 

.1450-. 1462 

.1702-. 1721 

.2151-. 2184 

.3258-. 3344 

19 

.1266-. 1281 

.1443-. 1462 

.1699-. 1721 

.2154-. 2184 

.3292-. 3359 

20 

.1271-. 1280 

.1451-. 1461 

.1706-. 1725 

.2162-. 2188 

.3296-. 3358 

30 

.1275-. 1286 

.1451-. 1470 

.1715-. 1733 

.2172-. 2197 

.3287-. 3355 

60 

.1276-. 1290 

.1459-. 1475 

.1727-. 1742 

.2200-. 2219 

.3310-. 3357 

100 

.1284-. 1297 

.1462-. 1479 

.1720-. 1742 

.2182-. 2212 

.3332-. 3392 

IS 


r 


TABLE  3 

MOMENTS  OF  THE  DISTRIBUTION  OF 

n 


M 

Mean 

Variance 

2 

.1124 

.00082 

.000026 

3 

,0929 

.00179 

.000107 

4 

.0898 

.00267 

.000281 

5 

.0892 

.00319 

.000409 

6 

.0894 

.00347 

.000471 

7 

.0899 

.00375 

.000550 

8 

.0897 

.00383 

,000565 

9 

.0903 

,00393 

.000550 

10 

.0902 

.00391 

.000546 

11 

,0906 

.00405 

.000599 

12 

.Ui/06 

.00405 

.000602 

13 

.0910 

.00410 

. 000609 

14 

.0908 

.00409 

.000603 

15 

.0909 

.00411 

.000619 

16 

.0911 

.00411 

.000599 

17 

.0914 

.00418 

.000643 

18 

.0915 

,00421 

.000639 

19 

.0914 

.00418 

.000617 

20 

.0914 

.00424 

.000625 

30 

.0917 

.00425 

.000632 

60 

.0920 

.00429 

.000629 

100 

.0922 

,00432 

. 000644 

«• 

.0926 

.00436 

.000640 

Figure  1.  Density  Function  of  C 
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